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ABSTRACT 

Apart from Mercury that has no known co-orbital companions, Venus remains as the inner 
planet that hosts the smallest number of known co-orbitals (two): (322756) 2001 CK 32 and 
2002 VE 68 . Both objects have absolute magnitudes 18 < H < 21 and were identified as Venus 
co-orbitals in 2004. Here, we analyse the orbit of the recently discovered asteroid 2012 XE133 
with H= 23.5 mag to conclude that it is a new Venus co-orbital currently following a transi- 
tional trajectory between Venus' Lagrangian points L 5 and L 3 . The object could have been a 
1 : 1 librator for several thousand years and it may leave the resonance with Venus within the 
next few hundred years, after a close encounter with the Earth. Our calculations show that its 
dynamical status as co-orbital, as well as that of the two previously known Venus co-orbitals, 
is controlled by the Earth-Moon system with Mercury playing a secondary role. The three 
temporary co-orbitals exhibit resonant (or near-resonant) behaviour with Mercury, Venus and 
the Earth and they follow rather chaotic but similar trajectories with e-folding times of the 
order of 100 yr. Out of the three co-orbitals, 2012 XEi 33 currently follows the most perturbed 
path. An actual collision with the Earth during the next 10000 yr cannot be discarded, an en- 
counter at 0.005 au may take place in 2028 but even closer encounters are possible within that 
time frame. Extrapolation of the number distribution of Venus co-orbitals as a function of the 
absolute magnitude suggests that dozens of objects similar to 2012 XEi 33 could be transient 
companions to Venus. Some additional objects that were or will be transient co-orbitals to 
Venus are also briefly discussed. 

Key words: celestial mechanics - minor planets, asteroids: general - minor planets, asteroids: 
individual: 2012 XE 133 - minor planets, asteroids: individual: (322756) 2001 CK 32 - minor 
planets, asteroids: individual: 2002 VE 6 8 - planets and satellites: individual: Venus. 



1 INTRODUCTION 

Early calculations carried out by Mikkola & Innanen (1992) 
showed that the existence of asteroids following co-orbital motion 
with Venus in the form of tadpole or horseshoe orbits stable for 
a few Myr was possible. Numerical integrations of the orbits of 
Near-Earth Objects (NEOs) revealed that some of them can even- 
tually follow co-orbital motion with Venus; for example, (4660) 
Nereus (Michel 1997) and (99907) 1989 VA (Christou 2000). In 
his study, Christou also predicted a quasi-steady-state flux of such 
objects moving in and out of the co-orbital regime. Tabachnik & 
Evans (2000) and Brasser & Lehto (2002) used simulations to 
show that Venus trojans would be stable in low-inclination orbits. 
These theoretical findings soon gained observational support with 
the discovery and subsequent identification of two present-day co- 
orbital companions to Venus: (322756) 2001 CK32, a horseshoe- 
quasi-satellite orbiter (Brasser et al. 2004), and 2002 VEes, a quasi- 
satellite (Mikkola et al. 2004). Both are relatively small and their re- 
spective absolute magnitudes, H, are 18.9 and 20.3. The existence 
of a transient population of NEOs trapped in a 1:1 mean motion 
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resonance with Venus has been explained by Morais & Morbidelli 
(2006) within the steady-state scenario envisioned by Bottke et al. 
(2000, 2002) where NEOs are constantly being supplied from the 
main asteroid belt. In their paper, it is predicted that the number of 
Venus co-orbital NEOs with H < 18 and H < 21 is 0.14±0.03 
and 1.6±0.3, respectively. Taking into account the two objects al- 
ready known, this translates into completeness for H < 21. If, as 
Morais & Morbidelli (2006) suggest in their work, current surveys 
have reached completeness at such small sizes, no objects in similar 
orbits and brighter than H = 21 should be observed in the future. 
On the other hand, numerical simulations indicate that the exis- 
tence of present-day primordial Venus co-orbitals is very unlikely 
due to the presence of multiple secular resonances (Scholl, Marzari 
& Tricarico 2005). 

In this paper, we show that the recently discovered asteroid 
2012 XE133 is a new transient Venus co-orbital with H = 23.5 
which is consistent with the Morais & Morbidelli (2006) predic- 
tions. The object was originally selected as a candidate because of 
its small relative semi-major axis, \a — avmus\ ~ 0.0004 au, and 
we use TV-body calculations to confirm its current co-orbital nature 
with Venus. This paper is organized as follows: in Section 2, we fo- 
cus on the available information on 2012 XE133, summarize some 
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Table 1. Heliocentric Keplerian orbital elements of asteroids 2012 XE133, (322756) 2001 CK32 and 2002 VEeg- Values include the lcr uncertainty. The orbit 
of 2012 XE133 is based on 102 observations with a data-arc span of 28 d. Although the quality of the orbit of 2012 XE133 is currently lower than that of the 
other two co-orbitals, it is comparable to the level of precision available when the other two objects were identified as orbital companions to Venus. (Epoch = 
JD2456400.5, 2013-Apr-18.0; J2000.0 ecliptic and equinox. Source: JPL Small-Body Database.) 



2012 XE 133 



(322756) 
2001 CK 32 



2002 VE 68 



Semi-major axis, a (au) 
Eccentricity, e 
Inclination, i (°) 
Longitude of the ascending node, Q (°) 
Argument of perihelion, u (°) 
Mean anomaly, M (°) 
Perihelion, q (au) 
Aphelion, Q (au) 
Absolute magnitude, H (mag) 



0.72297±0.00014 
0.4332±0.0003 
6.711 ±0.008 
281.096±0.006 
337.078±0.007 

30.98±0.06 
0.4098±0.0003 
1.0361±0.0002 
23.4±0.6 



0.725476965±0.000000004 
0.3824157±0.0000002 
8.13191±0.00002 
109.5025±0.0002 
234.1055±0.0002 
42.31827±0.00009 
0.44804316±0.00000011 
1 .0029 1 0769±0.0000()0006 
18.7±1.0 



0.7236506697±0.00000()0005 
0.4103603 1 ±0.00000005 
9.006018±0.000013 
231.582435±0.0000()5 
355.462382±0.000014 
163.09449±0.00003 
0.426693 16±0.00000()04 
1 .020608 1 798±0.00000()0007 
20.5±0.4 




Figure 1. The resonant angle, A r for the nominal orbit of 2012 XE133 in 
Table [T] (thick red line) and two of the control orbits. The thinnest blue 
line corresponds to a particular control orbit that was chosen close to the 
3cr limit so its orbital elements are most different from the nominal ones. 
The remaining orbit (green line) has osculating elements within \a of those 
of the nominal orbit. The output time-step for these plots is 1.1 d and no 
filtering or smoothing was applied to the data (same for subsequent figures 
unless explicitly noted). 



facts on co-orbital minor bodies and briefly outline our numerical 
model. Section 3 presents an analysis of the past, present and future 
dynamical evolution of 2012 XE133. Section 4 provides a compar- 
ative analysis that puts the object into context among the other two 
previously known Venus co-orbitals. Our results are discussed in 
Section 5 emphasizing the possible existence of additional objects 
like 2012 XE133. Our conclusions are summarized in Section 6. 



2 INITIAL CONDITIONS AND NUMERICAL MODEL 

Minor planet 2012 XE133 was discovered by Jess A. Johnson work- 
ing for the Catalina Sky Survey at Mt. Bigelow, Arizona on 2012 
December 12. Its apparent magnitude was V = 18.5. The object 
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Figure 2. The motion of 2012 XE133 over the time range (-200, 50) yr is 
displayed projected onto the ecliptic plane in a coordinate system rotating 
with Venus. The orbit and the position of Venus are also indicated. 



was observed on the same night by the Lunar and Planetary Labo- 
ratory (LPL)/Spacewatch II project and from the Magdalena Ridge 
Observatory, Socorro (McMillan et al. 2012Q With a value of the 
semi-major axis a = 0.7230 au, very close to that of Venus (0.7233 
au), this Aten asteroid is an NEO moving in a quite eccentric, e 
= 0.43, and slightly inclined, i — 6?7, orbit that makes it a Mer- 
cury grazer, Venus crosser, and Earth crosser. Therefore, its orbit is 
similar to those of the two previously known Venus co-orbitals (see 
Table[T): (322756) 2001 CK 32 and 2002 VE 68 . Its current orbit is 
based on 102 observations with a data-arc span of 28 d. In com- 
parison, the data-arc spans of (322756) 2001 CK32 and 2002 VEes 
were 512 d (with 82 observations) and 24 d (with 203 observa- 
tions), respectively, when they were identified as Venus co-orbitals 
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by Brasser et al. (2004) and Mikkola et al. (2004). As a recent dis- 
covery, little is known about this object besides its orbit with the 
exception of its absolute magnitude, H = 23.4±0.6 mag (assumed 
G = 0.15), which indicates a diameter in the range 62-138 m for 
an assumed albedo in the range 0.20-0.04. One of the objectives 
of this work is to encourage further observations of this small but 
interesting object. 

We selected 2012 XE133 as a candidate Venus co-orbital be- 
cause of its small relative semi-major axis, \a — av mus \ ~ 0.0004 
au. In the Solar system, an object is co-orbital to a planet if their or- 
bital periods and therefore their semi-major axes are nearly equal. 
Once an object is trapped in a 1:1 mean motion resonance and de- 
pending on its relative energy with respect to the host planet (Henon 
1969), it can describe any of the three main orbit types: quasi- 
satellite, tadpole or horseshoe. The names refer to the actual shapes 
of the orbits when seen in a frame of reference rotating with the 
host planet. The key parameter to differentiate among these three 
main orbit types and additional compound states, as first described 
by Namouni (1999), is the libration of the difference between the 
mean longitudes of the object and its host planet or relative mean 
longitude, A r . The mean longitude of an object is given by A = 
M + Q + uj, where M is the mean anomaly, Q is the longitude 
of ascending node and uj is the argument of perihelion. If the ob- 
ject follows a quasi-satellite orbit, A r oscillates around 0° . When 
the object moves in a tadpole orbit, A r librates around 60° (L4 tro- 
jan) or -60° (or 300°, L5 trojan). For objects moving in eccentric 
and/or inclined orbits the location of the Lagrangian points L4 and 
L5 changes; for example, in an orbit as eccentric as that of 2012 
XE133, it is shifted in A towards the L3 point by 30°-40° (Namouni 
& Murray 2000). If the libration amplitude is larger than 180°, en- 
compassing L3, L4 and L5 but not reaching the actual planet, it is 
said that the object follows a symmetric horseshoe orbit. These are 
the three principal dynamical states included in general co-orbital 
motion (see, e.g., Murray & Dermott 1999 and Mikkola et al. 2006 
for additional details). Compound states are also possible in which 
the object may librate around 0° with an amplitude > 180° en- 
compassing L4 and L5 (compound quasi-satellite-tadpole orbit), 
asymmetric horseshoe orbits (horseshoe-quasi-satellite orbiters) in 
which the libration amplitude > 270°, encompassing the planet, 
and a few other combinations (see, e.g., Namouni 1999; Namouni, 
Christou & Murray 1999). When the relative mean longitude oscil- 
lates freely, we say that the object is no longer in a 1 : 1 mean motion 
resonance with the planet, i.e. it becomes a passing object. 

In order to confirm the co-orbital nature of 2012 XE133 with 
Venus, we have performed TV-body calculations using the Her- 
mite integrator (Makino 1991; Aarseth 2003) in both directions of 
time for 10 kyr. Our model Solar system includes the perturbations 
by the eight major planets and treat the Earth and the Moon as 
two separate objects. It also includes the barycentre of the dwarf 
planet Pluto-Charon system and the three largest asteroids. The 
standard version of this direct TV-body code is publicly available 
from the IoA web sit^] Our integrations of the full equations of 
motion neglect relativistic terms as point (constant) mass objects 
moving in a conservative system are assumed. The role of the 
Yarkovsky and Yarkovsky-O'Keefe-Radzievskii-Paddack (YORP) 
effects (see, e.g., Bottke et al. 2006) is also ignored. We use ini- 
tial conditions (positions and velocities in the barycentre of the 
Solar system referred to the JD2456400.5 epoch) provided by the 
JPL's HORIZONS ephemeris system (Giorgini et al. 1996; Stan- 



dish 1998). In all the figures, t = coincides with the JD2456400.5 
epoch. Relative errors in the total energy at the end of the simula- 
tions are < 1 x 10~ lj . In addition to the calculations completed 
using the nominal orbital elements in Table[T|for 2012 XE133, we 
have performed 50 control simulations using sets of orbital ele- 
ments obtained from the nominal ones within the accepted uncer- 
tainties. The computed set of control orbits follows a normal distri- 
bution in the six-dimensional space of orbital elements. These syn- 
thetic orbital elements are compatible with the nominal ones within 
the 3a deviations (see Table[TJ and reflect the observational uncer- 
tainties in astrometry. The results of these calculations clearly show 
(see Section 3 and Fig. [TJ that the true phase-space trajectory will 
diverge exponentially from that obtained from the nominal orbital 
elements in Table [T] within a very short time-scale. This object and 
the other two Venus co-orbitals have e-folding times of the order of 
100 yr. The source of the Heliocentric Keplerian osculating orbital 
elements and uncertainties is the JPL Small-Body Databas^] Ad- 
ditional details can be found in de la Fuente Marcos & de la Fuente 
Marcos (2012). 



3 ORBITAL EVOLUTION: PAST, PRESENT AND 
FUTURE 

The motion of 2012 XE133 over the time range (-200, 50) yr as seen 
in a coordinate system rotating with Venus projected onto the eclip- 
tic plane is plotted in Fig. [2] and the three-dimensional evolution of 
its orbit in multiple frames of reference over the next few decades 
is displayed as an animation in Fig. [3] The object is a Venus co- 
orbital currently following a transitional trajectory between Venus' 
Lagrangian points L5 and L3. Our calculations (see Figs [T] and [4] 
left-hand panels) show that 2012 XE133 has already remained in the 
1:1 commensurability with Venus for several thousand years. For 
most of that time (4000 yr), it has been following a tadpole orbit 
in the neighbourhood of Venus' L5 point. After a close encounter 
with the Earth at 0.012 au 346 yr ago, it started moving away from 
the Lagrangian point L5 to undergo a transition between its previ- 
ous tadpole orbit and a horseshoe orbit. Horseshoe orbits are gen- 
erally less stable than tadpole orbits, and they are often considered 
as a transitional state before ejection from a 1:1 resonance (see, 
e.g., Murray & Dermott 1999). However, for the particular case of 
Earth's and Venus' horseshoe librators, long-lived configurations 
are possible (Cuk, Hamilton & Holman 2012). All the computed 
control orbits confirm that 2012 XE133 currently follows a rather 
irregular path of the horseshoe type, approaching the neighbour- 
hood of L3. The current precession rate of the bean-shaped loops 
associated with the asteroid orbital period is very slow and this tran- 
sitional path will last nearly 150-300 yr. Then the object's resonant 
angle will circulate after leaving the 1 : 1 resonance with Venus as a 
result of a very close encounter with the Earth at about 0.003 au in 
2098 CE. Its future orbital evolution depends strongly on the out- 
come of a very close encounter with the Earth that will take place in 
the time frame 85-340 yr from now. Nearly half of the control orbits 
evolve into passing orbits beyond that point, to never return to the 
Venus co-orbital region. The rest remain in or near the co-orbital 
region for thousands of years, at least ~3000 yr. For those orbits, a 
variety of dynamical behaviours are observed, including compound 
quasi-satellite-tadpole episodes and L5 tadpole orbits but also short 
circulation events. As for the past, nearly 70 per cent of the control 



2 http://www.ast.cam.ac.uk/~sverre/web/pages/nbody.htm 
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Figure 3. Three-dimensional evolution of the orbit of 2012 XE133 for a few decades starting on JD2456400.5, 2013-Apr-18.0, in three different frames of 
reference: heliocentric (left), frame corotating with Venus but centred on the Sun (top-right), and centred on Venus (bottom-right). The grey point denotes 
2012 XE133, the pink one denotes Mercury, the khaki one denotes Venus, the turquoise one denotes the Earth and the Sun appears in yellow. The osculating 
orbits are outlined and the viewing angle changes slowly to facilitate visualizing the overall orbital evolution. The current dynamical state resembles a tadpole 
orbit even if 2012 XE133 already left the Lagrangian point L5 after a close encounter with the Earth nearly 350 yr ago. 



orbits remained as an L5 trojan for about 3-5 kyr prior to the close 
encounter with the Earth that sends the object in the neighbourhood 
of L3. The remaining control orbits exhibit various horseshoe orbit 
types including not only symmetric but also asymmetric ones. Be- 
yond about 3000 yr into the past, the orbit is difficult to predict. 
Nearly 10 per cent of the control orbits remain within the Venus 
co-orbital region for the entire simulated time, 20 kyr, but alternat- 
ing multiple dynamical states and circulation episodes. Repetitive 
episodes as the ones described here, in which the relative mean 
longitude librates for several cycles and then circulates for a few 
more cycles before restarting libration once again, are character- 
istic of a type of dynamical behaviour known as resonance angle 
nodding, see Ketchum, Adams & Bloch (2013). These authors con- 
cluded that nodding often occurs when a small body is in an exter- 
nal (near) mean motion resonance with a larger planet. This type 
of complicated dynamics has been observed in other horseshoe li- 
brators. During the studied timespan, 2012 XE133 remains at a safe 
distance from Venus (see Fig.[5|, but close encounters are possible 
both at perihelion (with Mercury) and at aphelion (with the Earth); 
these encounters are critical in changing the dynamical status of the 
object as discussed in the following section. 



4 ASTEROID 2012 XE133 IN CONTEXT 

Fig.|4]displays the comparative evolution of the osculating orbital 
elements of all the known Venus co-orbitals (nominal orbits in Ta- 
ble [TJ. The short term three-dimensional evolution of the orbit of 
(322756) 2001 CK32 and 2002 VEes is displayed as an animation 
in Fig. [6] There are many similarities but also obvious differences, 
mostly in terms of short-term stability. All of them move in rather 
eccentric orbits with moderate inclinations and all of them are Mer- 
cury grazers, Venus crossers and Earth crossers. As the orbits of 
322756, 2002 VEsg and 2012 XE133 are very much alike, the over- 
all details of the mechanism that controls their orbital behaviour 
must be similar too. When moving in a 1:1 mean motion reso- 
nance, transfers between quasi-satellite, horseshoe and tadpole or- 
bits are the result of the libration of the nodes (Wiegert, Innanen & 
Mikkola 1998). Our analysis of the dynamics of 2002 VE68 in de la 
Fuente Marcos & de la Fuente Marcos (2012) confirmed the result 
obtained by Mikkola et al. (2004): during the quasi-satellite phase 
the distance to the descending node of the object remains remark- 
ably close to the value of Earth's aphelion and the distance to the 
ascending node is also relatively close to Mercury's aphelion. Here 
we have a similar situation. In the Solar system and for an object 
following an inclined path, close encounters with major planets are 
only possible in the vicinity of the nodes. The distance between the 
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Figure 4. Comparative time evolution of various parameters for the three known Venus co-orbitals: 2012 XE133 (left), (322756) 2001 CK32 (center) and 2002 
VEeg (right). The distance from the Earth (panel A); the value of the Hill sphere radius of the Earth, 0.0098 au, is displayed. The resonant angle, A r (panel 
B) for the nominal orbit in Table[T] The orbital elements a — ay cnus (panel C), e (panel D), i (panel E) and u> (panel F). The distance to the descending 
(thick line) and ascending nodes (dotted line) appears in panel G. Earth's and Mercury's aphelion and perihelion distances are also shown. The distance to the 
descending node of 2012 XE133 currently coincides with Earth's perihelion, that of (322756) 2001 CK32 is close to Mercury's aphelion and 2002 VEgs's 
matches Earth's aphelion. The output time-step for these plots is 0.3 yr. 



Sun and the nodes is given by r = a(l — e 2 )/(l ± e cosu), where 
the '+' sign is for the ascending node and the '-' sign is for the de- 
scending node. Close encounters at the nodes control the dynamics 
of 2002 VE68 and Fig. [4] panels A, B and G, clearly reveals that it 
is also the case for 322756 and 2012 XE133. In particular, the fu- 
ture evolution of the distances to the nodes of 2002 VEgg and 2012 
XE133 is analogous. We can therefore conclude that 322756, 2002 
VEe8 and the recently discovered 2012 XE133 are part of a dynam- 
ical group which is characterized by the libration of their nodes 
between the orbits of Mercury and the Earth. As a result, the grav- 
itational perturbations from the Earth and Mercury are currently 
most effective in keeping these asteroids at a safe distance from 



Venus. At present, 2012 XE133 approaches the Earth in the vicinity 
of its descending node in November-December every eighth year; 
in a way, this contributes to stabilizing the orbit if the distance of 
the closest approach is larger than about 0.03 au. These three ob- 
jects orbit the Sun in a near 8:13 resonance with the Earth but they 
also move in a near 9:23 resonance with Mercury although this is 
less clear than the near resonance with the Earth. The three objects 
exhibit resonant (or near resonant) behaviour with Mercury, Venus 
and the Earth. However, 2012 XE133 follows the most eccentric 
orbit and has the smallest inclination, which translates into more 
frequent and closer encounters with the Earth and explains why its 
orbit is the most unstable of the three. Encounters with Mercury as 
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Figure 5. The distance of 2012 XE133 from Mercury (top panel), Venus 
(middle panel), and the Earth (bottom panel). Close encounters are possible 
both at perihelion (with Mercury) and at aphelion (with the Earth). The 
values of the Hill sphere radii are also displayed. The output time-step for 
these plots is 1.1 d. 



close as 0.001 1 au and with Venus (only in the future) at 0.002 au 
have also been observed in a few of the control orbits. Close en- 
counters with Venus are only observed when the nodes are far from 
the perihelion/aphelion positions. 

Our results for 322756 are fully compatible with those in 
Brasser et al. (2004, fig. 1) for the next 2000 yr but differences 
appear beyond that point. These arise from the fact that we are us- 
ing the newest orbital elements, different physical models and the 
intrinsic chaos that drives the dynamical evolution of the object. 



5 DISCUSSION 

Morais & Morbidelli (2006) studied the population of Venus co- 
orbitals that are also NEOs and concluded that the average duration 
of a Venus co-orbital episode is ~32 000 yr. Our calculations in- 
dicate that 2012 XE133 current co-orbital episode is going to be 
significantly shorter but it is also true that its absolute magnitude is 
higher than the limiting value in their study and, therefore, their 
results may not apply in this case. They also estimated that the 
number of objects with absolute magnitude < H that become co- 
orbital with Venus over a time interval At is NAt/Atia, where 
N is given in table 1 of their 2006 paper and Ati-i is the aver- 
age duration of a Venus co-orbital episode. The model in Morais 

6 Morbidelli (2006) is based on a previous model by Bottke et al. 
(2000, 2002) which only applies in the size range corresponding to 
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Figure 7. The number distribution of Venus co-orbitals as a function of 
the absolute magnitude. This plot displays results obtained by Morais & 
Morbidelli (2006) with error bars, the current value of the counts (including 
2012 XE133) and the predictions from equation 



H < 22. This limit (about 400 m in diameter) reflects the transition 
from tensile strength dominated bodies to rubble piles or gravity- 
dominated aggregates (see, e.g., Love & Ahrens 1996). Even with 
this limitation in mind, it could be useful to have some estimate of 
the expected number of smaller objects. If we assume a power law 
of the form N(H) oc 10 aH and use the results obtained by Morais 
& Morbidelli (2006) to find the constant and the index a, we obtain 



N(H < Ho) = 6.98 x 10" 8 10° 



(1) 



where N(H < H a ) is the number of objects with absolute mag- 
nitude < H . This simple extrapolation (plotted in Fig. [7]) predicts 
about 14 objects for H < 24, which translates into about 10 objects 
similar to 2012 XE133 waiting to be discovered. But if the average 
duration of co-orbital episodes is shorter for smaller, less massive 
objects as the case of 2012 XE133 appears to suggest, then the ac- 
tual number of very small objects in co-orbital motion with Venus 
over a given timespan could be much larger in relative terms than 
the equivalent figure for larger objects like 322756 or 2002 VE68- It 
may be argued that because our present calculations neglect the ef- 
fects of non-gravitational forces like the Yarkovsky force, we can- 
not discuss the possibility of size-dependent trends in our results 
but nominal orbits used to compute the initial conditions for our 
integrations are based on the real positions of the objects, affected 
or not by non-gravitational forces. Therefore, the present orbital 
state of any observed object is always the result of all the (real) 
forces acting on the object. If the Yarkovsky force has any effects 
on the orbital evolution of the objects discussed here, those effects 
are already present on the initial conditions of the simulations but 
they are gradually lost if the integrations do not account for the 
Yarkovsky force, which is the case here. The time-scale in which 
such information is lost is likely longer than the characteristic e- 
folding time for these objects. 

The co-orbital identified in this paper was initially selected be- 
cause of its small relative semi-major axis that is nearly 0.0004 au. 
The relative semi-major axes of 322756 and 2002 VEes are 0.0021 
au and 0.00032 au, respectively. During the last two decades, wide- 
field NEO surveys have found thousands of objects passing in the 
neighbourhood of the Earth: could it be possible that some of these 
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Figure 6. Short term three-dimensional evolution of the orbit of (322756) 2001 CK32 (top) and 2002 VEs8 (bottom) in three different frames of reference: 
heliocentric (left), frame corotating with Venus but centred on the Sun (top-right) and centred on Venus (bottom-right). The purple point indicates 2012 XE133, 
the lime point indicates 2002 VEgg, the pink one indicates Mercury, the khaki one indicates Venus, the turquoise one indicates the Earth and the Sun appears 
in yellow. The osculating orbits are outlined and the viewing angle changes slowly to facilitate visualizing the overall orbital evolution. Asteroid 322756 is a 
horseshoe-quasi-satellite orbiter that currently follows the quasi-satellite stage of the compound orbit. 2002 VE68 exhibits typical quasi-satellite behaviour. 



objects may have not been properly identified and they are, in 
fact, trapped (even if temporarily) in a 1 : 1 commensurability with 
Venus? In order to explore this possibility and try to uncover ne- 
glected Venus co-orbitals, we have studied the evolution of all the 
known objects with relative semi-major axis (to Venus) <0.006 au. 
In this way, we have confirmed that no other known object (as of 
2013 March 7) is currently co-orbital to Venus. However, we have 



found two additional objects, 2007 AG and 2002 LT24, that may 
experience very brief co-orbital episodes (duration <10 3 yr) in the 
timespan ±10 4 yr. The asteroid 2007 AG has an absolute mag- 
nitude of 20.13 (it is brighter than 2002 VE6s), a relative semi- 
major axis of 0.0030 au, an eccentricity of 0.37 and an inclina- 
tion of 1 1?94. Its tadpole (L5) and horseshoe episodes are few and 
brief (< 900 yr); its orbit is well known with a data- arc span of 
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Figure 8. Same as Fig.[4]but for (99907) 1989 VA. The nodes match Mer- 
cury's perihelion periodically and for extended periods of time. This minor 
body likely was a Venus co-orbital in the past. 

1791 d. More interesting is the case of 2002 LT 24 at H = 21.85 
with a relative semi-major axis of 0.0035 au, an eccentricity of 
0.5 and an inclination of just 0?76. It also experiences horseshoe 
behaviour with respect to Venus in several occasions, the longest 
centred about 2000 yr into the future and lasting nearly 1000 yr. 
The orbit of this object is also very reliable with a data-arc span of 
2934 d. Like 2007 AG, 2002 LT 24 is an Aten and NEO that is also 
included in the Potentially Hazardous Asteroid (PHA) list. Finally, 
2003 KO2 appears to signal the edge of the co-rotational region of 
Venus at a relative semi-major axis of 0.0041 au. Again, an Aten, 
NEO and PHA with H = 20.1 1, an eccentricity of 0.5 and an incli- 
nation of 23? 5, it experiences a single half horseshoe loop during 
the studied timespan. All the objects move in rather eccentric or- 
bits. In summary and although the model in Morais & Morbidelli 
(2006) merely provides statistical expectations, currently available 
data are compatible with predictions in Fig. [7] 

Outside the timespan ±10 4 yr, we confirm the early result ob- 
tained by Namouni et al. (1999) and Christou (2000) and pointed 
out above, that the asteroid 99907 (relative semi-major axis 0.0051 
au) was a Venus co-orbital. Christou (2000) found that in the ±10 
yr interval this object has a 50 per cent probability of following 
co-orbital motion with Venus; we find that the object had multi- 
ple co-orbital episodes in the past, some of them lasting thousands 
of years. However, its dynamics (see Fig. [8} is somewhat differ- 
ent from that of the objects previously discussed as it moves in a 
more eccentric (0.59) and inclined path (28? 8). Now the distance 
of one of the nodal points of 99907 from the Sun usually coincides 
with the minimum distance to Mercury from the Sun. Every 40 
kyr both nodes reach Mercury's neighbourhood; this translates into 
more frequent close fiy-bys with the innermost planet around that 
time (see Fig. [9}. 

Apart from Mercury that has no known co-orbital compan- 
ions, Venus still remains as the inner planet that hosts the smallest 
known number of co-orbitals. From an observational point of view, 
the scarcity of detected co-orbitals for the two innermost planets 
is not surprising and it does not necessarily translate into a real, 
physical lack of objects. If the majority of the actual innermost co- 
orbitals have H > 21 as equation {T} suggests and spend most of 
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Figure 9. Same as Fig.|5]but for 1989 VA. 

the time far from the Earth and away from the ecliptic plane fol- 
lowing orbits similar to those of the three known Venus co-orbitals, 
their apparent magnitude may be too faint to be detected by the av- 
erage telescope used to discover NEOs even at the most favourable 
geometry; the maximum elongation when the object's angular sep- 
aration from the Sun is the largest and it is best observed from 
Earth. It is not by chance that 2012 XE133 was discovered in 2012 
December as it was its closest approach to the Earth since 1940 
December when it passed 0.05 au from our planet. It had a close 
encounter with the Earth on 2012 December 19 at 0.055 au and it 
passed 0.048 au from Mercury on 2013 March 23. On the other 
hand, if most Venus co-orbitals are as faint as 2012 XE133 but they 
do not move in paths that get close to Earth's (for example, low ec- 
centricity co-orbitals), their detectability by ground-based surveys 
is severely compromised. Only space-borne programmes may be 
successful in this case. 

Asteroid 2012 XE133 is part of the list compiled by the JPL 
Sentry Syste nfjthat includes objects that may be involved in future 
Earth impact events. Our calculations indicate that an actual colli- 
sion with the Earth during the next 10 000 yr cannot be completely 
ruled out. On 2028 December 30, 2012 XE133 may pass close to 
the Earth and the Moon at a distance < 0.006 au from both bod- 
ies. A similar encounter will take place in 2049; even closer will be 
the encounter in 2098 at 0.003 au. In 2151 April, it will approach 
0.0027 au from Earth. Some of our control orbits record close ap- 
proaches below 0.00068 au. 



4 http://neo.jpl.nasa.gov/risk/ 
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6 CONCLUSIONS 

In this paper we have identified a new high eccentricity Venus co- 
orbital, 2012 XE133. This identification is consistent with Morais 
& Morbidelli (2006) predictions in terms of absolute magnitude. 
Our analysis of the orbit of 2012 XE133 reveals that this small 
NEO is following a transitional, highly chaotic path. The object has 
been a co-orbital companion to Venus for several thousand years 
and it may become a passing object within the next few hundred 
years. The future orbital evolution of this recently discovered mi- 
nor body is not easily predictable beyond that time-scale. Its orbit is 
the most chaotic among Venus co-orbitals. The three known Venus 
co-orbitals follow similar orbits with e-folding times of the order 
of 100 yr and exhibit resonant (or near-resonant) behaviour with 
Mercury, Venus and the Earth. Our calculations indicate that the 
thickness of Venus' co-orbital region is 0.004 au with 2007 AG, 
2002 LT24 and 2003 KO2 signalling the edge of the region. Extrap- 
olation of the number distribution of Venus co-orbitals as a function 
of the absolute magnitude suggests that dozens of objects similar to 
2012 XE133 could be transient companions to Venus but ground- 
based surveys may have substantial difficulties in finding them due 
to the nature of their orbits and intrinsic faintness. 

Our calculations did not include the Yarkovsky effect which 
may have a non-negligible role on the medium, long-term evolu- 
tion of objects as small as 2012 XE133. Proper modeling of the 
Yarkovsky force requires knowledge on the physical properties of 
the objects involved (for example, rotation rate, albedo, bulk den- 
sity, surface conductivity, emissivity) which is not the case for 2012 
XE133. Detailed observations during future encounters with the 
Earth should be able to provide that information. The non-inclusion 
of this effect has no impact on the assessment of the current status 
of this object as a Venus co-orbital. 



Christou A. A., 2000, Icarus, 144, 1 

Cuk M., Hamilton D. P., Holman M. J., 2012, MNRAS, 426, 305 1 

de la Fuente Marcos C, de la Fuente Marcos R., 2012, MNRAS, 427, 728 

Giorgini J. D. et al., 1996, BAAS, 28, 1 158 

Henon M., 1969, A&A, 1, 223 

Ketchum J. A., Adams F. C, Bloch A. M., 2013, ApJ, 762, 71 

Love S. G., Ahrens T. J., 1996, Icarus, 124, 141 

Makino J., 1991, ApJ, 369, 200 

McMillan R. S. et al., 2012, MPEC 2012-X85 

Michel P., 1997, A&A, 328, L5 

Mikkola S., Innanen K., 1992, AJ, 104, 1641 

Mikkola S., Brasser R., Wiegert P., Innanen K., 2004, MNRAS, 351, L63 
Mikkola S., Innanen K., Wiegert P., Connors M., Brasser R., 2006, MN- 
RAS, 369, 15 

Morais M. H. M., Morbidelli A., 2006, Icarus, 185, 29 

Murray C. D., Dermott S. F., 1999, Solar System Dynamics, Cambridge 

Univ. Press, Cambridge, p. 97 
Namouni F, 1999, Icarus, 137, 293 

Namouni F, Murray C. D., 2000, Celest. Mech. Dyn. Astron., 76, 131 
Namouni F, Christou A. A., Murray C. D., 1999, Phys. Rev. Lett., 83, 2506 
Scholl H., Marzari F., Tricarico P., 2005, AJ, 130, 2912 
Standish E. M., 1998, JPL Planetary and Lunar Ephemerides, 

DE405/LE405. Interoffice Memo. 312.F-98-048, NASA JPL 
Tabachnik S., Evans N. W., 2000, MNRAS, 319, 63 
Wiegert P. A., Innanen K. A., Mikkola S., 1998, AJ, 1 15, 2604 



ACKNOWLEDGMENTS 

The authors would like to thank S. J. Aarseth for providing the 
code used in this research and the referee for his/her constructive 
and useful reports. This work was partially supported by the Span- 
ish 'Comunidad de Madrid' under grant CAM S2009/ESP-1496 
(Dinamica Estelar y Sistemas Planetarios). The authors also thank 
M. J. Fernandez-Figueroa, M. Rego Fernandez and the Department 
of Astrophysics of the Universidad Complutense de Madrid (UCM) 
for providing excellent computing facilities. Most of the calcula- 
tions and part of the data analysis were completed on the 'Servidor 
Central de Calculo' of the UCM and the authors thank S. Cano 
Alsiia for his help during this stage. In preparation of this paper, 
the authors made use of the NASA Astrophysics Data System, the 
ASTRO-PH e-print server and the MPC data server. 



REFERENCES 

Aarseth S. J., 2003, Gravitational N-body simulations. Cambridge Univ. 

Press, Cambridge, p. 27 
Bottke W. F, Jedicke R., Morbidelli A., Petit J. M., Gladman B., 2000, Sci, 

288, 2190 

Bottke W. F, Morbidelli A., Jedicke R., Petit J. M., Levison H., Michel P., 

Metcalfe T., 2002, Icarus, 156, 399 
Bottke W. F., Jr., Vokrouhlicky D., Rubincam D. P., Nesvorny D., 2006, 

Annu. Rev. Earth Planet. Sci., 34, 157 
Brasser R., Lehto H. J., 2002, MNRAS, 334, 241 

Brasser R., Innanen K. A., Connors M., Veillet C, Wiegert P., Mikkola S., 
Chodas P. W, 2004, Icarus, 171, 102 



